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■ Abstract 
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Upon passing an a.c. electrical current along magnetic micro- or nanostrips, the measurement 
of a d.c. voltage that depends sensitively on current frequency and applied field has been recently 
reported by A. Yamaguchi and coworkers. It was attributed to the excitation of spin waves by 
the spin transfer torque, leading to a time-varying anisotropic magnetoresistance and, by mixing 

o 

O . of a.c. current and resistance, to a d.c. voltage. We have performed a quantitative analysis 

i \ by micromagnetics, including the spin transfer torque terms considered usually, of this situation. 

The signals found from the spin transfer torque effect are several orders of magnitude below the 

00 : 

I/*} ■ experimental values, even if a static inhomogeneity of magnetization (the so-called ripple) is taken 

into account. On the other hand, the presence of a small non-zero average (Ersted field is shown 

qq ' to be consistent with the full set of experimental results, both qualitatively and quantitatively. 

O ■ 

We examine, quantitatively, several sources for this average field and point to the contacts to the 

^ • 

^ . sample as a likely origin. 

PACS numbers: 72.25.Pn, 76.50.+g, 41.20.-q, 72.15.-v 
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I. INTRODUCTION 



The possibility to act on the magnetization of a sample by an electrical current within it, 
not through the classical (Ersted field but through the spin-polarization of electrical current 
in ferro magnets, offers fascinating opportunities in nano magnetism and nanoelectronics [l, 
Q . In the situation where the sample consists of separated and uniformly magnetized media 
crossed by the current, the description of the physics appears simpler and, indeed, agreement 
between experiments and modelling does not appear out of reach {3, 4]. However, when the 
current flows in a magnetic medium with a continuously varying magnetization, the situation 
is more complex. As a result, several forms for this so-called spin transfer torque (STT) 

8], and the appropriate equation for magnetization dynamics 

m 

In such a situation, the more experimental results in different configurations is clearly 
the better. Among these, the recent discovery of an electrical rectification effect in magnetic 
strips with widths of the order of a micrometer and thicknesses of the order of a few tens of 
nanometers is especially appealing. The effect was observed for current densities below 
or of the order of those required for STT to act on domain walls. However, a relatively 
large static field was applied so that the strip was in a single domain state, contrarily to 



have been proposed 0, O, [7 
has even been questionned [9 



the situation where a signal was measured in presence of a domain wall [12|, [13]. This 
last feature is puzzling. Indeed, STT within a continuous magnetization structure is only 
expected when a magnetization gradient exists. In the simplest STT formulation, valid for 
slow magnetization variations with respect to electrons' spin precession or diffusion length, 
the STT is namely expressed as 

dm 
~dt 



\STT 



u ■ V ) rh + (3fh x 



u ■ V I rh 



(1) 



where the velocity u is an expression of the current density J with spin polarization P 
according to 



u 



J 



9VbP 
2eM ' 



(2) 



The (small) number f3 has been related to spin flip of the conduction electrons, in several 



models 



0, y, q, q. 



From (Cp), one sees that magnetization gradients along the electric 
field are required. Such gradients should however not exist in the experimental situation 
considered above (long strip under a large field), at least for perfect samples. The possibility 
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FIG. 1: Geometry of the sample and notations definition. The sample geometry with notations 
definition is shown in (a). The two calculation models are schematically depicted: (b) periodic 
model and (c) infinite model. The hatched area depicts the current x profile. 

mentioned by the authors is that the uniform state becomes unstable under a.c. current at 
an appropriate frequency, as indeed predicted for very large d.c. currents flrj| |. 

The object of this paper is to perform a full micromagnetic analysis of the situation in 
order to analyze the various sources of rectification signal discussed above, and to quantita- 
tively compare the calculated si gnal s with the experimental results. 

The experimental conditions |ll|] are as follows (see Fig. [1] for notations): the sample is 
a magnetic strip, several micrometers long with various widths (from 300 to 5000 nm) and 
thicknesses (30 to 50 nm), with the experimental constraint of a close to 50 Q resistance; a 
magnetic field H app is applied in the sample plane, at an angle #h; an a.c. current with swept 
frequency is injected into the sample through a coplanar waveguide; current densities are of 
;he order of 10 10 A/m 2 i.e. low compared to those required for domain wall displacement 



171 ] . The experimental results may be summarized by: (i) a d.c. voltage is measured with a 
marked frequency dependence, it becomes important (~ fiV) only at a well defined frequency 
of the order of the ferromagnetic resonance frequency; (ii) the position of the resonance is 
strongly influenced by the field magnitude (in accord with Kittel's law); (Hi) the d.c. voltage 
increases as the square of the injected current; (iv) the angle dependence of the d.c. voltage 
is well described by a sin(26 l H) cos(6 , H ) law at large fields, turning to sm(9u) at low fields. 

Our approach uses both analytical and numerical micromagnetics: we solve the Landau- 
Lifchitz-Gilbert magnetization dynamics equation with incorporation of the two basic STT 
terms (0Q). Section [11] is devoted to the case of a uniform magnetization along the strip 
axis, where no STT is expected in the linear limit. The next section introduces structural 



inhomogeneities, in the spirit of the well known ripple patterns in thin films 18], so that 
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a STT is present in the ground state. As both configurations lead to d.c. voltages much 
below the experimental levels, Sec. [IV] investigates the effects of a small average (Ersted 
field. Finally, as the samples have been carefully designed to avoid fields from the current 
leads, an intrinsic origin to this field from different electron scattering properties at the top 
and bottom surfaces of the strip is discussed. 



II. UNIFORM MAGNETIZATION 



We first look at the simplest situation where the magnetization does not change in the 
direction of the electric field, at least in the rest state without current. 



A. Analytical analysis of the uniform situation 

The magnetization at rest will be denoted m , with \rho\ = 1. In the presence of the a.c. 
current, a small deviation rh(f,t) appears (\rh\ <C 1). The current is described by a spatially 
uniform u that is harmonic in time with pulsation to. With the axes defined in Fig. [T]one has 
u = (u(t) = u cos(u;t), 0, 0). The LLG equation supplemented by both STT terms reads, to 
first oder in the deviations m 



dm ( -± -> _ \ _ dm 

— = 7o yH xm + hx m J + am x — 



dirt dirt 

— u Q cos(ut)— h Buq cos(u)t)mQ x — — . (3) 

ox ox 

In this equation, Hq is the effective field of the static magnetization (with Hq x m = 
by definition) and h is the effective field resulting from the existence of the deviation m, 
with contributions from the exchange and magnetostatic energies. With only the first two 
terms on the right-hand side of ([3]), upon diagonalization, the various spin wave modes 
corresponding to the static magnetization m are obtained jT l| . Their amplitude is fixed 
by the thermal noise. The last two terms represent the spin-wave pumping by the STT. 
We note that they have the form of a product u(t)m (we forget the spatial derivatives here 
as the argument is about the time dependence). Therefore, if fh varies in time with the 
pulsation uj also, these terms do not contribute to m as they lead to pulsations 2uo and 0. 

In fact, a coupling of the form u(t)m is known to give rise to parametric excitation, 
i.e. the generation of a solution at frequency / by pumping at frequency 2f. Parametric 
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pumping of spin waves from a uniform starting configuration would not be in agreement 
with the experimental results (point (i)), since the resonant frequency found for the current 
is of the order of the FMR frequency, not twice this value. 



B. Numerical calculations 

As analytical calculations suffer from some limitations, such as linearization and the 
consideration of simple structures only, they were completed by full micromagnetic numerical 
simulations (see Methods in Ref. j^oj]). These were performed by solving the LLG equation 
with the STT terms. The typical velocity uq equivalent to the current applied was uq = 
3.25 m/s (corresponding to the experimental current density J = 6.5 x 10 10 A/m 2 with a 
polarization P = 0.7). In addition, for the purpose of showing better the effect of STT, 
values of uq as large as 100 or even 1000 m/s were applied. As no effects of the (3 term were 
observed (this will become clear from the analytical calculations), the results shown below 
were obtained with (3 = 0. We will throughout the paper show results for just one sample 
size (width w = 300 nm and thickness t = 50 nm), i.e. values corresponding to one sample of 
Ref. l"j. Material parameters were those representative of NigoFe2o, namely magnetization 
M s = 800 kA/m, exchange energy constant A = 1 x 10~ n J/m, no crystalline anisotropy 
(except for the 'ripple' case, see Sec. IHIj) . gyromagnetic ratio 70 = 2.21 x 10 5 m/(As) 
and damping parameter a = 0.01. A static field H &vv> was applied in the sample plane, 
at an angle 9h, with standard values poH app = 40 mT and 9h = 45° (the value where 
the d.c. signal is close to maximum). The calculation region length L, a part of the real 
sample length L s , was taken to be L = 1 u,m, mostly (calculations with L = 2 or 4 fim 
were also conducted for the purpose of checking the dependence on L of the d.c. voltage). 
The mesh size was 4 x 4 x 50 nm 3 mostly. The d.c. voltage was computed from the 
time variation of the anisotropic magnetoresistance (AMR) of the sample. Denoting the 
resistivity change upon magnetization rotation by Ap (with, for the used NiFe alloy at 
room temperature, Ap ^ 0.5 x 10~ 8 Qm), the dependence of the wire resistance upon its 
magnetization distribution is expressed as 

AMR = -^<ml + m 2 z >, (4) 
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where S is the wire cross-section area (S = wt) and <> denotes the average over the 
calculation region. 

Note that, as there is no domain wall in the calculation region and a field with a transverse 
component is applied, the calculation scheme has to be different from that used for the 
simulation of domain wall dynamics 20J: one cannot assume that outside the calculation 
region the magnetization is uniform and equal to (±1,0,0). Thus, the calculation region 
was embedded in a wire of infinite length according to two different models (too simple 
embedding schemes can lead to gradients of fa in the x direction, that cause large spurious 
spin transfer torques). In the 'periodic' model (Fig. Hb), the calculation region is supposed 
to be repeated periodically in the x direction. The exchange and demagnetizing fields 
are calculated accordingly, as well as the x gradients for the STT. The current density is 
uniform. In the 'infinite' model (Fig. [lb), one assumes that the values at the left edge of 
the calculation region extend to infinity on the left side, and similarly on the right. The 
boundary conditions at left and right are free. The demagnetizing field takes into account 
these two semi-inifinite regions. In order to avoid end effects, the current density is zero at 
x = 0, L and rises (along a length L/3) linearly towards the set value at the center of the 
calculation region. 

As seen in the analytical analysis (Sec. Ill Ap . no effect is expected in first order per- 
turbation as the initial state magnetization has no gradient along the current direction. 
Numerically however, the avoidance of any x gradient is impossible, as this would require an 
infinitely precise numerical evaluation of the demagnetizing and exchange fields. It follows 
that, depending on the model used and the accuracy of the numerical scheme, a gradient 
along x of the magnetization remains that gives rise to a non-zero STT, therefore to a 
magnetization oscillation and finally to some d.c. voltage. The most uniform initial state 
was obtained with the infinite model, and is depicted in Fig. [2] for a field H &pp applied at 
9h = 45°. Note that, as the field is inclined and its value /io-f/app = 40 mT is below the effec- 
tive transverse anisotropy of the magnetic strip, the magnetization becomes non uniform in 
the transverse y direction, with more rotation towards the field at the strip y center. Fig. [2] 
proves that the non uniformity in the x direction of the magnetization is essentially due to 
the numerical precision of the calculations (a double precision number is described with a 
relative precision of 10~ 16 ). 

The time variation of the magnetization at position (x = L/2, y = w/A) in the calculation 
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FIG. 2: Maps of the initial magnetization state under a static field /io-^app = 40 mT {H &pv = 
400 Oe) applied in the film plane at an angle Oh = 45°, in a perfect wire 300 nm wide and 50 nm 
thick. The axial component m x is shown in (a), with the gray scale applied to the deviation of m x 
from its average, magnified 50 times. For the m y component (b), the gray scale is magnified 20 
times (from 0.1 to 0.2). In order to display the magnetization (non) uniformity in x despite the non 
uniformity in y due to the inclined applied field, the maps of the x differential (m(i+l, j) — m(i— 1, j), 
with i and j cell indices) of the magnetization components are shown in (c) for m x and (d) for m y , 
magnified by a factor 10 15 . The images cover the calculation region length (1 //m). 

region, under a.c. current at various frequencies, is plotted in Fig. [3] for the initial state 
depicted in Fig. [2J The current amplitude was u$ = 100 m/s, much higher than the largest 
experimental value (~ 3 m/s), but this was necessary in order to see some effect. Indeed, 
the deviations of this central moment are extremely small, of the order of 10~ 15 whereas the 
initial value is m y = 0.151 at that point. The Fourier analysis of the data reveals that the 
fundamental frequency fo of the a.c. current is seen in the spectra only when it is close to 
a resonant frequency of the sample (the numerical calculation of the FMR-like excitation 
of the sample by an a.c. field in this configuration shows resonance at /o — 11 GHz). In 
this regime, multiples of the fundamental frequency are also seen. In addition, when / is 
close to twice this resonant frequency, a (broader) response at half the current frequency can 
be seen. This is a characteristic of the parametric excitation, expected from the analytical 
analysis presented above. 

In order to see how this oscillation behaves in space, some snapshots of the magnetization 
structure are provided in Fig. HI The images reveal that the complex wavy state seen in 
the structure at rest (Fig. [2]) acquires a tiny oscillation in the form of a quasi- standing wave 
pattern (with some very geometrical features that may be artifacts). Such a pattern had 



been evidenced in the preliminary calculations [l_l|, [2l|. However, the amplitude of this 
oscillation is extremely small (in cases where the x gradients of static magnetization at the 
left and right edges of the calculation region are less carefully avoided, the oscillations are 
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FIG. 3: (color online) Time variation of the magnetization at one cell in the calculation region, 
under a.c. currents of varying frequency and fixed amplitude uq = 100 m/s, in a perfect sample 
(to numerical accuracy). The transverse component m y is shown, for the sake of sensitivity. Real 
time traces are shown in (a), once the steady state oscillation has been reached, for two frequencies 
/o of the a.c. current (note the 10 15 magnification). In (b) Fourier spectra are displayed, for 
lower frequencies (fo = 5 to 15 GHz with 1 GHz step, top panel) and higher frequencies (/o = 19 
to 24 GHz with 1 GHz step, bottom panel; note the magnified vertical scale). They reveal the 
excitation at the fundamental frequency with, in some frequency range, the apparition of under- 
multiple peaks when they correspond to characteristic frequencies of the sample. The current 
frequencies /o corresponding to some spectra are indicated. 

stronger) . 

Therefore, the computed d.c. voltages are, even at their maximum around the FMR 
frequency, orders of magnitude (here, ~ 15) below the experimental values, as shown in 
Fig. [5j Thus, clearly, another mechanism has to be found. 



III. NON UNIFORM MAGNETIZATION 

The analysis in the preceding section has shown that one of the reasons for the absence of 
signal was that the initial state contained no gradient along the electric field. A well-known 

a result of a random 



cause for such non uniformitiy is the so-called 'ripple' pattern [18| 
distribution of anisotropy in crystallites. The response of the magnetization structure to 
this random potential is to organize itself in ripple domains, where the magnetization de- 
viates slightly from its average value in one direction or the other. These domains have a 
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FIG. 4: Images of the magnetization oscillating state under a.c. current at /q = 11 GHz, with 
uq = 100 m/s, in a perfect sample. The first image is built from the centered x differential of the 
magnetization y component, magnified by a factor 10 15 , at time t = 4 ns after current application 
(in the steady regime, note however the closeness to Fig. [2ti). The other images are generated 
by difference of the magnetization y component with the same component at time t = 4.0 ns, 
magnified by 3 x 10 14 . 
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FIG. 5: Spectrum of the d.c. voltage for the uniform initial state. The current amplitudes are 
very large, no = 100 and 1000 m/s. The values for no = 1000 m/s are divided by a factor 10 2 in 
order to compare them to those at no = 100 m/s. Note the extreme smallness of the d.c. voltage 
obtained (vertical scale in atto- volts (10 -18 V). 



characteristic lens shape, with the long axis of the lens oriented normal to the magnetiza- 
tion. The presence of the ripple structure is attested by Lorentz electron microscopy (see 



e.g. Ref. 



m). 



A. Analytical model for smooth ripple 

A first model consists in assuming that the ripple pattern can be described by a smooth 
undulation of the magnetization. Thus we write mo(x,y) = (cos(6 l o), sin(6 l o), 0), where 9q 
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oscillates, more or less regularly, in space. The small deviation m is decomposed into an 
in-plane component and an out-of-plane component: m = ae$ + be z , where eg is the unit 
in-plane vector orthogonal to the local magnetization and e z the unit vector normal to the 
film plane. The linearized LLG equation projected on these two vectors now reads 

da db , .d9 

Tt = loHb ~ a di~ u{t) ^' 

KAJ V KAlV kj Ju 

db rT da „ , . d9n . . 

~dt ~ ~dt ~dx (5) 

To solve the equation, we assume that the two components of the effective field are simply 

proportional to the magnetization deviations : H a = —h a a and Hb = —h^b. For a thin fim, 

one expects that hb « M s if the ripple period is much larger than the sample thickness. 

The term h a corresponds to the effective potential that stabilizes the value of 6q and its 

variation with x. The equations are linear and contain the STT as a driving term now. 

They are easily solved in the harmonic approximation : a = Aexp(iut), b = B exp(zu;t) and 

u = U exp(icjt). We then get 

A-U— — + - ^° hb + aiUJ ^ (6) 

dx (1 + a 2 )uo 2 — ^Ihahb — aiu(h a + h b ) 

The magnetization in-plane deviation becomes large at the ferromagnetic resonance defined 

by uj 2 cs (1 + a 2 ) = 7o^a^-fe- As the meaning of the deviation component a is a rotation of the 

in-plane magnetization angle Q , we see that at resonance a current u cos(o; res t) transforms 

the structure as 

u cos(u ms t) 



0o(x) -> 



x — 



(7) 



aio{h a + h b ) 

This means that the ripple magnetic pattern under STT is set into oscillation along the x di- 
rection. The amplitude of oscillation at resonance is typically (in the thin film approximation 
where h a <§C hb ~ M B ) 



Ax w — (8) 

For uq = 3.25 m/s, the oscillation amplitude computed from (|Sj) is Ax = ±2 nm. 

The predicted oscillation of the ripple pattern is therefore quite small. Moreover, in an 
infinite wire, the modulation of the total anisotropic magnetoresistance is exactly zero, as 
the structure is merely translated. For a wire of finite length, some small signal can be 
expected as, at both ends, some part of the structure disappears or appears because of the 
structure oscillation in position. In such situation however, the sign of the d.c. voltage 
should be arbitrary and the d.c. voltage magnitude should not depend on sample length. 
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FIG. 6: The initial magnetization state for a sample with ripple structure. The magnetization x 
component is displayed in (a) with, similarly to Fig. [2] a magnification by a factor 50. In order to 
better display the magnetization variations in x despite the non uniformity in y due to the applied 
field, the map of the x differential of m y is shown in (b), where the centered cell to cell difference 
was magnified (slightly over) by a factor I0 4 . The m y (x) profile averaged over ten lines of cells 
at the strip center is shown in (c), with its Fourier transform in order to display the longitudinal 
'periodicity' of the ripple pattern. 

B. Numerical calculations 



As for the preceding section, numerical calculations were performed in order to get quan- 
titative results. For inducing a ripple structure, a random anisotropy field was introduced, 
with a fixed value and a random in-plane easy axis orientation. Results obtained for 
^oHk — 3 mT only will be shown here. The description of the magnetic structure at rest is 
provided in Fig. [61 to be compared to that of the perfect sample (Fig. [2]). It is clear that 
the random anisotropy leads to much larger magnetization deviations than the residual ones 
in the perfect case. The structure evidences some periodicity along x (about 250 nm, see a 
profile of m y in the x direction together with its Fourier transform in Fig. Eb). In the ripple 
theory 18] , two structural periods appear, that in the direction of magnetization being con- 
sequently smaller than that in the orthogonal direction. This larger period appears to be 
mostly suppressed by the reduced width of the nanostrip sample, for the parameters chosen 
here. Note the similarity of the differential magnetization image in Fig. [6)d with typical 
TEM images in Lorentz mode, that also depend on magnetization gradients 221 ] . 

Under application of an a.c. spin polarized current, the magnetization structure is driven 
into some oscillation. Fig. [7] provides some snapshots of the magnetization distribution, for 
the high current amplitude u = 100 m/s. Despite the large current, the oscillation nature is 
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hard to apprehend from the movies of m y (t). Therefore, the time differentials m y (t) — m y (to) 
were drawn. The y component of fh was chosen as it is the most labile, the magnetization 
at rest under the applied field being close to the x direction, in the standard case illustrated 
here. A very smooth modulation appears, that contrasts with the noisy appearance of the 
rest state gradient image. It is a standing wave pattern tuned to the current frequency (this 
pattern is already discernable in Fig. H] behind the more geometrical texture). Note that the 
wave pattern is not exactly 'standing': the modulation breathes with time but also deforms 
slightly during one period. This deformation looks like the pattern oscillation computed in 
the previous subsection, but this requires a more precise verification. The 'wavevector' of 
this pattern is roughly parallel to the x direction, and the wavelength is decreasing as the 
applied frequency increases (not shown). This last feature, compared to the conclusions of 
the previous subsection, shows that the dominant mechanism of the spin transfer torque 
action is different. 

A more appropriate description can be constructed by assuming that the ripple pattern 
is random. The spin transfer torque due to the x gradients of the structure at rest m is 
equivalent to a field that 'pumps' the deviations m, expressed as 



Note that Fig. [6b can also be seen as displaying the largest component of this field. Therefore, 
in first approximation, this field is random in space but perfectly harmonic in time. By 
Fourier transformation, this field can excite spin wave modes with a frequency matching the 



is similar to the thermal field that gives a non zero amplitude to the thermodynamic spin 
waves (as seen in Brillouin scattering), with the only difference that instead of being white 
in temporal frequency it is monochromatic. Therefore, the observation of induced spin wave 
modes is no surprise. From the analogy with the thermal field we get that the amplitude of 
the spin wave is proportional to current (and also to the strength of the ripple drriQ y /dx). 

The next question concerns the detection of this induced spin wave as a d.c. voltage. 
In first approximation, the wavelength of the standing spin wave pattern (A ~ 60 nm 
at 11 GHz here) being much smaller than the length of the sample, the variation of the 
total AMR averages to zero. However, as the standing wave pattern is not perfect, being 
deformed by the ripple structure (see Fig. [7]), the average is not perfectly zero. It is a random 




(9) 



current frequency f$ and a wavevector contained in the spectrum of H r . Thus, the field H r 
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FIG. 7: Images of the magnetization oscillating state under a.c. current at frequency /q = 11 GHz, 
with no = 100 m/s, for a sample with ripple structure. The first image, at time t = 8 ns, is the 
centered x differential of the magnetization y component, magnified by a factor 10 3 . It is very 
close to the structure at rest (less magnified). The next images, corresponding to one period of 
oscillation, are just the differences of m y (t) with m y (t = 8 ns), also magnified by a factor 10 3 . 

number with a typical value X/L times smaller than the local amplitude. Fig. [8] shows a 
spectrum calculated for one particular realization of the ripple. A resonance appears (with 
a close to expected shape for that spectrum) at /o — H GHz, i.e. the FMR frequency of 
the sample, meaning that the uniform mode is also pumped by H r . The values are much 
larger than in the nominally uniform case studied earlier (Fig. ED and a large irregularity 
is present, a signature of the randomness of the averaging over the sample size. To this 
curve is superposed the average of the absolute values of the d.c. voltages calculated for 16 
realizations of the ripple. Furthermore, a calculation run for a L = 2 fim calculation length 
and u = 100 m/s gives a similar result, as expected (not shown). Finally, the voltage is found 
to be proportional to the square of the a.c. current, as expected from the proportionality 
of spin wave amplitude to current (compare the scaled spectra in Fig. [H]for u = 100 and 
1000 m/s). This has the consequence that, for current amplitudes similar to the experiments, 
the calculated d.c. voltage should be smaller by a factor 1000, i.e. a million times smaller 
than what is measured. Therefore, we conclude that the ripple hypothesis is also not able 
to account quantitatively for the measured rectification signal. 



IV. EFFECT OF A PARASITIC FIELD 



The analysis of the first sections has shown that virtually no rectification effect is to be 
obtained through a STT as described in (CQ), in a single domain ferromagnetic wire. Of 
course, we cannot exclude that other mechanisms than those already investigated play a 
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FIG. 8: Spectrum of the d.c. voltage for one initial state with ripple. The current amplitude is 
very large, no = 100, and the length of the calculation box is L = 1 um. The spectrum obtained 
with uq = 1000 m/s, divided by a factor 100, is superposed (dashed curve). The dash-dot curve is 
an average, over 16 realization of a ripple pattern, of the absolute signals (their sign being random). 

role in the interpretation of the experiments, but we need also to investigate the effect of 
mechanisms competing with STT. In the following, we examine the influence of a parasitic 
(Ersted field on the rectification effect. 

If we consider an infinitely long wire with a uniform current density, the average value of 
the (Ersted field is zero. As the wire has a flat cross-section, the largest components of the 
field are vertical (z direction) and are found at the two lateral edges of the wire. For the 
300 x 50 nm 2 wire submitted to a current density J = 6.5 x 10 10 A/m 2 , the maximum vertical 
field is equal to n$H z = ±2.3 mT (the horizontal components reach /j,qH v = ±1.8 mT). These 
fields, normal to the magnetization at rest, excite magnetization oscillations in the vicinity 
of the wire surfaces. What will be the impact on the sample resistance ? As a static field 
is applied in the sample plane, at angle 9h from the wire axis, the m y component at rest is 
not zero. In presence of the (Ersted field, this component oscillates and therefore the AMR 
is modulated. However, as the fields have opposite signs at two opposite surfaces, the total 
AMR modulation should cancel by integration. 

Note that, in principle, the cancellation may be incomplete. For example, in the spin- 
wave community, the effect of interfaces is often taken into account by a pinning length £ 
such that the boundary condition for the oscillating magnetization fa is 

dm m 

± ftT + I = ' (10) 



where n denotes the coordinate normal to the interface 
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24|. The original condition 23] 



specified that £ = A/K s , with K s the surface anisotropy constant, but recently magnetization 
oscillation profile calculations were shown to be explainable through a pinning length related 
to the sample dimensions 24j . actually a pure magnetostatic effect. Here, whereas there is 
no reason to have any dissymetry between opposite surfaces for magnetostatics, a chemical 
or structural difference between the top and bottom surfaces of the magnetic film cannot 
be ruled out, leading to different surface anisotropies. It is obvious that this will result in a 
non zero AMR oscillation and, therefore, a non zero d.c. voltage at the resonance frequency 
of this n = 1 perpendicular standing spin wave mode (PSSW). Generally, the frequency of 
the PSSW is different from that of the uniform FMR mode, allowing their discrimination. 
Indeed, on the one hand, frequency should increase due to exchange. But, on the other 
hand, the transformation of the lateral magnetostatic dynamic charges from monopolar to 
dipolar decreases the dynamic demagnetizing field, thus reducing the frequency. In the case 
at hand, the direct numerical calculation gave 10.6 GHz for the PSSW, close to the value of 
11 GHz for the uniform mode, so that a frequency discrimination is not possible. Turning 
now to the signal levels, an upper bound to this contribution can be obtained by evaluating 
the d.c. voltage on one half of the sample thickness. As peak values of ±1.25 /iV were 
obtained, we conclude that the d.c. voltage due to a perfect CErsted field is not sufficient, 
by a factor ~ 10, to explain the experimental results Let us therefore now suppose that the 
average of the CErsted field is not exactly zero. 



A. Analytical model with a.c. field excitation 

We assume here that the magnetization is uniform (in-plane angle 9q) so that from (jl]) the 
total AMR is -ApL s sm 2 (6 )/S. When the magnetization angle is oscillating as 8q + a {t), 
with \a(t)\ << 1, the sample resistance change due to AMR, for a sample of length L s , reads 

R(t) « -ARsm{26 )a{t). (11) 

where AR = ApL s / S. We now determine the angle oscillation a(t). The equations of motion 
use the same variables as in Sec. IIII Al but now instead of considering the STT we include 
an a.c. field h = (0, h y , h z ) (this represents a general non-zero average for the CErsted field 
due to the current, and neglects the zero average part). The LLG equations of motion now 
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read 



da db 

— = -j h b b - a — + j h z , 
dt dt 

db , da 

— = 7 /i a a + a— - j hy cos(6> ). (12) 

The harmonic solution (with h y = H y exp(iu)t) etc.) reads similarly 

A _ iuioHz + (loh + aiu)joH y cos(6> ) 

u 2 (l + a 2 ) - "Yohhb - aioj^oiha + h b ) ' 

At resonance, taking into account that h b 3> h a , the amplitude is approximately (note that 
the y field gives the largest effect) 

^_.7oiWflo). (14) 

auj vcs 

As this variation of magnetization does not change sign over the sample, it will not average 
to zero by integration. Consequently, this magnetization oscillation produces an oscillation 
of the sample resistance at frequency u, from which a d.c. voltage results. The d.c. voltage 
has a peak close to the resonance frequency : as the phase of A is close to 7r/2 the signal 
is zero at resonance but has peaks of alternating signs on both sides (as in the experiments 



ll|). The typical d.c. voltage is thus 



^ yP = ^sin(2^)cos(^)^. (15) 

With this formula we see that, as in the experiments, the voltage is proportional to the 
current square (as the (Ersted field is proportional to the current), and that the dependence 
on magnetization angle conforms to the sin(26 l ) cos(6 l ) law. Fig. [9] displays the frequency 
dependent d.c. voltage evaluated from ffTTl [T5|) . The d.c. voltages are of the order of the 
micro- Volt, like in the experiments. With a sole a.c. field in the y direction, an antisymmetric 
line shape is obtained, whereas an a.c. z field contributes with a symmetric line shape. 



Ill ] may be interpreted by a combination of both 



Therefore, the experimental line shapes 
of these fields. 

We conclude analytically that a non-zero (Ersted field provides a plausible explanation 
of the experimental results showing a rectification effect in nano and microstrips. This was 
qualitatively stated in Ref. {2]]], but dismissed on the ground that the average field should 
be zero. The last section of this paper will therefore study possible origins for such fields. 
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FIG. 9: Plot of the frequency dependent d.c. voltage according to (fl~3|) and (fTT|) . The frequency is 
normalized to resonance, hf, = M s and h a is such that resonance occurs at 11 GHz, a = 0.01, and 
the magnitude of the average a.c. field is HoH y = 0.1 mT. The dash and dot lines were computed 
by adding an a.c. field along the z direction equal to ± the y field. This changes the peak shape 
(less antisymmetric) but not the signal p.p. amplitude, to first order. The voltage ARIq/2 is 
163 fiV for Jo = 6.5 x 10 10 A/m 2 in the 300 x 50 nm 2 nanostrip of length 1 /zm, so that the d.c. 
signal is calculated to be ±1.4 (j,V. 

B. Numerical calculations 

The influence of the (Ersted field was investigated by numerical calculations also. Indeed, 
the magnetization non uniformity in the transverse direction, as a function of applied field 
angle and magnitude, or the evaluation of the effect of the full (Ersted field, require a 
micromagnetic simulation. 

For the latter point, calculations with a variation along the not essential. Thus, 

for the description of the full (Ersted field in the sample cross-section and the calculation 
of the resulting d.c. voltage, a 2D model (y, z) with invariance in the x direction was also 
employed. The case with no bias (for h y or h z ) was already discussed. When a bias is added, 
we find that the d.c. voltage spectra are affected by the presence of the full (Ersted field only 
at the frequency of the PSSW, with a minor quantitative influence when the d.c. voltage is 
computed on the full sample thickness. Therefore, we kept the model used in the rest of the 
paper for the other calculations with a.c. field, and neglected the full (Ersted field whose 
average is zero. 

First, an a.c. field with y component only was considered h y = H y cos(ut) with an 
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FIG. 10: Computed d.c. voltage spectra obtained for increasing values of the applied field, for two 
field angles 45° and 225°, and sample dimensions 300 x 50 x 1000 nm 3 . An a.c. field of amplitude 
fioHy = 0.1 mT is taken into account as an average (Ersted field, and the voltage is computed for 
an a.c. current density Jq = 6.5 x 10 10 A/m 2 . The static field values are chosen to be those of the 



experiment 



lil. 



amplitude HoH y = 0.1 mT. The dependence of the d.c. voltage spectra on the static field 
value, for two opposite field directions, is illustrated in Fig. [TQ1 The field values and angles, as 
well as the current density, were chosen to match those of the experiments [11] . In comparison 
to the analytical calculation assuming uniform magnetization (Fig. [H]), the signal is roughly 
divided by 2. The results are very similar to experiments, qualitatively and quantitatively 
(the computations apply to a calculation region 1 /im long whereas for these dimensions and 
the nominal resistivity of NiFe the sample length corresponding to 50 Q is 3 /im), with only 
a difference in the peak shape. However, from the analytical modeling (see Fig. EJ) we know 
that such change of shape can be obtained by adding a z component to the a.c. field. 

Another important experimental feature is the dependence of the d.c. voltage peak to 
peak amplitude on the applied field angle. It was analytically shown above that the observed 
behaviour (the sin(2^n) cos(#h) law) was expected for an a.c. h y field. The numerical results 
substantiate this conclusion, as shown in Fig. [TT] for two values of the static applied field. 
For the large value (/io-^app = 0.5 T), the magnetization angle 6q follows the field angle #h 
closely and the law is well obeyed. For the smaller field (fioH app = 0.16 T), closer to the 
shape anisotropy field of the nanostrip, an evolution towards a sin^n) variation is evident. 

Therefore, assuming that a non zero average (Ersted field exists, we have been able to 
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FIG. 11: Computed d.c. voltage maximum amplitude as a function of d.c. field angle. In order 
to saturate the nanostrip, a large field has to be applied (/io-^app = 0.5 T here). A lower field 
(/Uo-f^app = 0.16 T, open squares) leads to a deformed curve. Other conditions are identical to 
those for Fig. [TU1 The solid line draws the cos(2## ) sin(#// ) law. 

reproduce all experimental results, quantitatively with a.c. fields ^iqH « 0.1 mT. This is in 
sharp contrast with the alternative explanations based on STT, that result in rectification 
signals that are orders of magnitude smaller. 



V. POSSIBLE ORIGINS FOR A PARASITIC FIELD 



The question now arises naturally about the origin of such a non zero average field. We 



note that the design of the electrical connections to the magnetic wire in the experiments 11] 
was very symmetrical, the sample being inserted in a coplanar waveguide with a ground- 
signal-ground structure. This ensures geometrically that the (Ersted field applied to the 
sample in the presence of a current is minimized. Of course, any imbalance in current 
backflow between the 2 ground leads, or imperfect centering of the sample, creates an a.c. z 
field (for example, we estimate roughly that a 1 % imbalance would give \iqH z ~ 0.1 /iT ). 
However, there are two different reasons for which a non-zero average y component of this 
field can still exist: (i) the gold coplanar waveguide (100 nm thick) contacts the magnetic 
sample by its top surface; (ii) due to differences in the top and bottom interfaces of the 
magnetic film, the current distribution could be non uniform, even in an infinitely long 
magnetic strip. We try now to estimate these fields. 
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A. Field from current distribution in the sample thickness 



We first assume that neither a non magnetic conductive underlayer nor overlayer exist 
(that would obviously give rise to a y field) , and look for an intrinsic origin for a non uniform 
current distribution. For the description of electronic transport in thin metallic films, the 



Fuchs-Sondheimer model [25j evaluates the current distribution in the thickness of a thin 
film, within a Boltzmann equation approach. A key parameter introduced in this model is 
the specularity parameter p at every interface: p = 1 means that electron reflection at the 
interface is perfectly specular and p = that it is random. In the latter case the current is 
reduced at the interface. The characteristic scale (in z) over which this reduction extends 
is given by the mean free path A. Extending the Fuchs-Sondheimer calculation to different 
interfaces at the film top (z = t) and bottom (z = 0) [26], we obtain in the limit of a thick 
film (t ^> A) the following current distribution 



J x (z) = a E x 



2 V A 7 2 A 
where the function $ is an exponential integral 



(16) 



$(x) = - [ exp(--)(l-u 2 )du. (17) 
2 Jo u 

Two z profiles of the current for the case t — 50 nm and a mean free path A = 2.5 nm are 
plotted in Fig. [121 The cases p = 0, p' = and p = 1/3, p' = 2/3 are compared. In the 
former, the current density is reduced to one half at both interfaces, symetrically. In the 
latter, the current is reduced to 2/3 at one interface and to 5/6 at the other, so that J x (z) 
becomes asymmetric. 

The next step is to evaluate, from this asymmetric current distribution, the average 
(Ersted field within the sample. As J x is asymmetric in z, this gives rise to a non zero 
average for the field component transverse to the strip axis (H y ). The analytical calculation 
of this average field shows that the contribution of a current layer at depth z to the average 
field is very close to linear, the contribution being zero for z = t/2. In other words, the field 
average over y and z is 

2 A r t/2 

< H y >= / J x (z')z'dz' (18) 

J -t/2 

where A is a number equal to n/2 in the limit w ^> t and amounting to A = 1.296 for the 
sample we consider here with w/t = 6. From ([TBI) and ([TBI we obtain for J = 6.5x 10 10 A/m 2 , 
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FIG. 12: Fuchs-Sondheimer model of the current distribution in a film with thickness t = 20A 
where A is the electron mean free path. The parameters p and p' are the specularity parameters 
for the electron reflection at the bottom and top surfaces, respectively. The two horizontal scales 
allow zooming at the vicinity of the top and bottom surfaces. 

w = 300 nm and t = 50 nm, ^H y « 5 fiT. This value could be perhaps doubled by taking 
a longer mean free path and increasingly different specularity parameters, but remains too 
low for explaining the experiments. 

B. Field from the contact regions 

In the lack of a full 3D current distribution calculation, the typical field due to the gold 
electrodes being deposited on top of the strip surface can be evaluated roughly. We assume 
that, at both ends of the magnetic strip (length L a ), the current flows vertically between 
the magnetic strip and the gold electrodes, the latter being twice as thick and 10 times 
more conductive than the magnetic sample. The length of this vertical part is taken to be 
the thickness t of the magnetic layer. At the center of the magnetic strip, considering that 
L s ^> t, the resulting typical y field reads 



With J = 6.5 x 10 10 A/m 2 , L s = 3 jum, w = 300 nm and t = 50 nm, one gets ^qH v ~ 5 fiT. 
Thus the central value of this field is also too low with regard to the experiments. However, 
it becomes much larger close to the contacts. For the average value of the field over the full 




(19) 
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sample (x and y), we obtain 

<H v >=^H2j) (20) 

(in order to remove a divergence we assume that the current flows vertically over a length t 
along the x axis at the contacts). The same numbers now give fio < H y >= 0.05 mT, a value 
compatible with those explaining the d.c. voltages experimentally measured. Therefore, we 
propose that the rectification signals are due to this field. It should be noted that this field 
decreases as the inverse of the sample length, in contrast to the intrinsic mechanism whose 
the contribution is independent of sample length. 



VI. CONCLUSION 



This work has tried to reach a quantitative description of the spectra of d.c. voltage 
versus frequency obtained on nanostrips subjected to a static in plane field. We have found 
that, for a perfect wire, no signal should be expected owing to a linear analysis. The 
consideration of internal magnetic inhomogeneities in the nanostrip (ripple structure) do 
lead to non zero rectification voltages. However, these are disorder dependent and much 
too low compared to experimental findings. As a solution to this paradox, we tested the 
influence of a non zero average value of the (Ersted field generated by the a.c. current, 
and found that small values (0.1 mT) of this average lead to quantitative agreement with 
experimental findings. We propose two origins for these fields. The extrinsic one is due to the 
contacts put on the magnetic sample (as considered recently in the case of vortex excitation 
by a.c. currents 27]]). The intrinsic one is due to an asymmetric current distribution in 
the thickness of the magnetic nanostrip. In the semi-classical description by Fuchs and 
Sondheimer of the conductivity of very thin films, such an asymmetry results from different 
specularity parameters for electron reflection at the top and bottom surfaces of the magnetic 
layer. 

Therefore, we conclude that (Ersted field effects should be carefully investigated when 
interpreting experiments by the spin transfer torque mechanism. The extrinsic contribution 
should be precisely evaluated, so that the numerical evaluation of the field taking into 
account the full sample structure may become mandantory in the future. The important 
role of small a.c. fields at resonant frequencies is not so surprising in fact: in ferromagnetic 
resonance, for a sample with a = 0.01 like NiFe, the excitation field is of the order of 
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HqH = 1 /iT. 
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[ I 

Note added in proof : A. Yamaguchi, in a recent work [28( , comes also to the conclusion 
that a field effect is very probably the origin of the rectification signal. In addition, as 
source of the non zero average field, he proposes the deformation of the field pattern at high 
frequencies due to electromagnetic effects, the permittivity of the substrate being different 
from that of air. 
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